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ABSTRACT 

KIC 1571511 is a 14d eclipsing binary (EB) in the Kepler dataset. The secondary of 
this EB is a very low mass star with a mass of 0.14136 ± O.OOOSGM© and a radius 
of O.17831^Q QQQg2i?0 (statistical errors only). The overall system parameters make 
KIC 1571511B an ideal "benchmark object": among the smallest, lightest and best- 
described stars known, smaller even than some known exoplanet. Currently available 
photometry encompasses only a small part of the total: future Kepler data releases 
promise to constrain many of the properties of KIC 157151 IB to unprecedented level. 
However, as in many spectroscopic single-lined systems, the current error budget is 
dominated by the modeling errors of the primary and not by the above statistical 
errors. We conclude that detecting the RV signal of the secondary component is crucial 
to achieving the full potential of this possible benchmark object for the study of low 
mass stars. 
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1 INTRODUCTION 

Current stellar models describe well the basic properties of 
a wide range of stars. However, there are persisting discrep- 
ancies at the lower end of the stellar mass range, where very 
late-type stars have measured radii that are higher than the- 
ory predicts (e.g. Lacy 1977, Ribas 2006, Torres & Ribas 
2002). This letter reports the discovery that Kepler target 
KIC 1571511 is an EB which contains a very low-mass star 
as secondary - hereafter just 51 IB - based on the public 
Q0-Q2 Kepler data and our radial velocity (RV) follow-up. 
We point out the special location of 51 IB in the parameter 
space and offer to use it as a benchmark object for future 
studies of low-mass stellar objects. In i|2] we highlight the 
important properties of the KIC 1571511 photometry and 
subsequent RV follow-up, in i|3]we model the system, in !|4] 



* Based on observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Fin- 
land, Iceland, Norway, and Sweden, in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica de 
Canarias 

f E-mail: avivofir@wise.tau.ac.il 



we estimate the system's physical parameters based on all 
available data, and conclude in i|5] 



2 IDENTIFICATION AND FOLLOW-UP 

We identified KIC 1571511 as an interesting system as early 
as four days after the Kepler QO-Ql data was made pub- 
lic [Borucki & the Kepler Team 2010, see references therein 
for a description of the Kepler satellite and data]. This sys- 
tem contains periodic eclipse-like events every 14.02d, with 
a depth of ~ 2% - a conspicuous signal at Kepler's high 
precision. At that time we suspected this system to host a 
giant transiting planet overlooked by the Kepler team due 
to its (assumed) orbital eccentricity. It is noteworthy that 
while high-eccentricity transiting planets can generate tran- 
sits that are "too" long and photometrically identical to EBs 
on a circular orbit, the Kepler pipeline does not consider ec- 
centric orbits in its Data Validation module El We thus were 

^ Kepler Data Processing Handbook §9.3, document number 
KSCL19081-001 of 1 April 2011. 
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Table 1. Measured Radial Velocities 
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eclipse variation, such as the reflection effect or the Doppler 
boosting signal (Zucker et al. 2007). We chose a strong fil- 
ter to better remove the stellar activity since it has higher 
amplitude than the the Doppler boosting signal and since 
anyhow no further information is expected from it given our 
high precision RV data. Once an initial solution was ob- 
tained (below) we rejected 22 points that deviated by more 
than 4(T from the model, and this data was used for the 
light-curve part of the final solution. 



able to secure a few RV measurements with the FIES spec- 
trograph at the NOT observatory to test the giant-eccentric- 
planet hypothesis. We note that this object was also consid- 
ered as an overlooked planetary candidate by Coughlin et 
al. (2010) but they did not provide RV data. 

RV follow-up of KIC 1571511 was started in October 
2010 using the FIES fiber-fed Echelle spectrograph attached 
to the 2.61 m Nordic Optical Telescope (observing program 
P41-426). The observations were performed under grey/dark 
time with good and stable sky-conditions. The 1.3 /im high- 
resolution fiber was employed, yielding a resolving power of 
A/AA ~ 67, 000 and a wavelength coverage of about 3600 - 
7400 A. Following the method described in Buchhave et al. 
(2010), long-exposed ThAr spectra were acquired immedi- 
ately before and after each target spectrum to improve the 
wavelength solution and trace any instrumental drifts. Stan- 
dard IRAF routines were used for the data reduction and 
spectra extraction. The RV measurements were derived by 
cross-correlating the target spectra with a spectrum of the 
RV standard star HD182488 (Udry et al. , 1999) observed 
with the same instrument set-up as the target. A journal of 
the FIES observations is given in Tablejl] The FIES observa- 
tions revealed a single-line spectroscopic binary (SBl) with 
an eccentric orbit (e = 0.33) and a RV semi-amplitude of 
K = 10.5 km/s, compatible with a very low-mass compan- 
ion star orbiting the main component. Figure [l] (right-hand 
panels) show the FIES RV measurements along with the 
Keplerian RV curve resulting from the best-fitting simulta- 
neous photometric-RV solution (see qSl and residuals. 



3 SYSTEM MODELING WITH MCMC 
3.1 Pre-processing 

Since both the primary and secondary eclipses were easily 
identified using a weak filter we were able to use a stronger 
filtering scheme for the final hght curve (hereafter LC) with- 
out modifying the eclipses themselves. We computed an iter- 
ative ~ 1 day long, second-order Savitzky-Golay filter with 
5a outlier rejection using only the out-of-eclipses sections of 
the data, and separately for each continuous section of the 
Q0-Q2 data. We then interpolated the filter to the times 
of the eclipses and normalized the LC with that filter. No 
data points were rejected to this point, but only removed 
from the filter calculation, and this LC contained 6177 data 
points. Continuous sections were taken between quarters 
and anomalies, as reported by the Kepler Release Notes 
(We considered the following anomalies: attitude tweak, safe 
mode. Earth point, and coarse point). We note that such 
a strong filter automatically removes any in-phase out of 



3.2 Methods 

We used MCMC (Monte Carlo Markov Chains, e.g. Tegmark 
et al. 2004] for the simultaneous solution of the photometry 
and RV data. As in many MCMC codes, an MCMC chain 
begins by computing the model for a given set of parameters. 
We then add a random perturbation to each parameter p 
chosen from normal distribution of width ap. If the total 
has reduced at the new perturbed location, the perturbed 
parameters set is accepted as the new set. If the total 
has increased, the new set is only sometimes accepted - at a 

2 _ 2 

probability of exp{ "^"'"^ ) ^ also commonly referred to as 
the Hastings-Metropolis jump condition. Early chains used 
a rather large ap, or step size, as the parameters space was 
explored for interesting regions of low x^- Once an initial 
solution was obtained we set all the different jump sizes 
to the standard deviation of that MCMC test-chain, and 

1/2 

re-scaled them down by a common factor of Nj^Qp where 
Ndof ~ 13 is the number of degrees of freedom. This is 
done since once all the parameters' jumps are measured in 
units of their own standard deviation, the MCMC process is 
actually a random walk in a TVooF-dimensional space, and 

1/2 

SO the typical distance covered by such a step is NjjQp. 

We use the above MCMC procedure and the Mandel 
and Agol (2002) formalism to model the system as two lumi- 
nous spheres (vs. Roche geometry models). We expand the 
formalism to account for the secondary's flux and Kepler's 
finite integration time (Kipping 2010). Specifically, eq. 40 
of Kipping (2010) implied a sub-sampling of A'' ~ 2.3 sub- 
samples to reduce the modeling errors to below the mea- 
surement errors. We therefore chose = 5 to make sure 
this effect is indeed minimized. 

Our model included eccentric Keplerian orbits with a 
period P, a unit semi major axis a = 1, and eccentricity and 
argument of periastron given by ecos(tj) and eain{ui). The 
fitted reference time parameter is the more easily (and more 
accurately) observed time of mid eclipse Tmid while the time 
of periastron passage is computed from it using the previ- 
ous parameters. Other parameters are the fractional radius 
of the primary r\/a, the relative radius of the secondary 
r2/r\ and the orbital inclination to the line of sight i. Us- 
ing the Mandel and Agol (2002) limb-darkening model we 
inluded a quadratic limb-darkening for the primary using 
■ui,i, «2,i and a linear limb-darkening for the secondary with 
■iti,2. Since the total flux is normalized to unity, we only var- 
ied the secondary fractional luminosity L2 / Ltot ■ When con- 
tamination, or "third light" L3, was included contamination 
levels of more than ~ 3% gave poorer flts, and contamination 
levels < 3% were just as good as the fit with zero contami- 
nation, so we adopted a fixed L3 = 0. Finally, the mass ratio 
Q = — -.'"^ affects the positions of the two bodies - and so 
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Table 2. Photospheric parameters of KIC 1571511A as derived 
from the spectral analysis of the co-added FIES spectrum 



Teff[° K] log(g) [cm s-2] Fe/H [dex] iisini [km s^^] 
6195 ±50 4.53 ±0.10 0.37 ± 0.08 7.9 ± 0.5 



it is included in the LC model too, and not just in the RV 
model. From the above it follows that for the RV model we 
only needed to add the systemic velocity 7 and the overall 
scale of the system - how many metres there are in one semi 
major axis a. These are only a shift- and scaling- terms for 
the otherwise known RV morphology (from the previous pa- 
rameters). In principle, the scaling parameter should be an 
MCMC variable and the systemic velocity should be fitted 
analytically - but since this is a single-lined binary the scale 
is degenerate with the mass ratio. We therefore did not vary 
the scaling parameter but fixed at a value computed from 
the estimated mass of the primary (Q, the above period 
and mass ratio, and Kepler's laws. 



3.3 Application to KIC 1571511 

We searched the above parameter space using numerous 
MCMC chains until we got very close (in retrospect - within 
Ax^ < 5) of the absolute minimum. We then ran 30 
chains of 5 x 10* steps each to densely sample the local vol- 
ume. We verified that all of them converged on the same 
parameters set and then concatenated all the chains (as 
they are independent, see Tegmark et al. 2004) to a final 
1.5 X 10® steps long chain used for the parameters value and 
errors estimation. We also checked that all parameters are 
well mixed, i.e., their effective length (which is their nominal 
length divided by their autocorrelation length) is ^ 1. Since 
we started all the runs from a point very near the final min- 
imum no "burn-in" was required and all steps were kept. To 
further check that no farther and deeper minima exist, we 
also ran 300 shorter chains that did not start from near the 
reported Xmin - but randomly perturbed by up to lOOa in 
every parameter. No other deeper minimum was found. The 
best-fitting parameters and their errors were determined as 
the median values of each parameter's distribution and the 
the ranges that span of the steps on either side of that 

median. These parameters and errors, as well as some de- 
rived quantities, are reported in Table [3] and the LC and RV 
models are shown on Figure [l] 

Overall the quality of the fit is satisfactory, with xLt of 
8046.5, or reduced x?ed of 1.31. Some of the excess residuals 
can be attributed to imperfect filtering of microactivity on 
primary: the EB is active with an amplitude of 7.2 x lO^** 
while the LC's residuals are 1.8 x 10~* - almost 4 times 
lower. It is unlikely that this variability is dominated by 
511B since this would imply a variability of about 25% in 
its flux on ~ day time scales, and so we attribute it almost 
entirely to the primary (hereafter just 511A). It is notewor- 
thy that while both limb-darkening coefficients of 511A were 
constrained by the data, no such constraint is yet possible 
even for the simplest linear model of 511B. 



4 PHYSICAL PARAMETERS ESTIMATION 

Single-lined spectroscopic and eclipsing binaries, such as 
KIC 1571511, do not allow for the full model-free deter- 
mination of their parameters. The derived quantities can 
only be solved up to a single line of possible mass/radius 
relations for each one of the components (Beatty et al. 
2007, hearafter B07). One therefore needs to derive the pri- 
mary's mass from some models, and systematic errors at this 
stage are the overwhelming source of error for this dataset 
(more below). We therefor took extra care and used mul- 
tiple tools redundantly for the following step. We modeled 
the co-added FIES spectrum of 511 A primarily using the 
new SPC fitting scheme (Buchhave et al. , in prep.) which 
allowed us to extract precise stellar parameters from the 
spectrum (see Table |2] We double checked this new anal- 
ysis with the more traditional spectral synthesis package 
SME (Valenti & Piskunov 1996; Valenti & Fischer 2005) 
and got a similar and consistent results. We then used Teg, 
log(g) and Fe/H and a grid of Yonsei-Yale model isochrones 
(Yi et al. 2001), and performed a Monte Carlo analysis 
to infer the stellar mass and radius and an estimate of 
their uncertainties. This yielded a staller mass and radius 
of M511A = 1.2651^30^^0 and 7?5iia = 1.216lg;J^^7?0 re- 
spectively. Since the radius of 511 A can also be inferred from 
the above mass and the obsrvationally-constrained relation 
of B07, and since the latter gives lower errors, we adopt 
its values: iisiiA = 1.3431q oio^q. We note that at this 
stage we effectively have three different determinations for 
the primary's surface gravity: one (log(g)=4.53) comes di- 
rectly from the spectral analysis, and the other two are in- 
direct from the combination of mass and radius from the 
above Monte Carlo distribution (log(g)=4.37) and B07 rela- 
tion (log(g)=4.28). We choose to adopt this last determina- 
tion because it is least model-dependent: given the primary's 
mass it only assumes Kepler's laws and spherical stars. 

511B is thus determined to have a mass of 0.14136 ± 
3.6 X 10~*; X lO-^Mo and a radius of 0.17831+6 2 x 
10"*; Ijg X lO~*i?0 (statistical and modeling errors, respec- 
tively). Not surprisingly, both the LC and the RV data are 
far more precise than the model-dependent derivation of the 
mass of 511 A, and the latter is the overwhelming source of 
error in the derived physical parameters of 511B. We there- 
fore quote two error estimations of the mass and radius of 
511B in Table[3]- the lar ger is derived only from the mod- 
eling error of the mass of 511A, and the smaller is derived 
only from statistical errors of the LC-RV fit. 

Since we did not model the reflection effect the 
L2/Ltot = 2.75 X 10~* depth of the secondary eclipse is a 
combination of the primary's flux reflected off the secondary, 
and of the intrinsic luminosity of 51 IB itself. The former 
can be calculated as Frefiected = ^sV^)"^ ^i'l^i}) where Ag 
is the geometric albedo and a is the primary-secondary 
distance during secondary eclipse, which in this case is: 
Fre fleeted = X 4.2 X 10~^. We find that refiected light can- 
not contribute more than about 1/6 of the light lost during 
the secondary eclipse. Thus at least 2.33 x 10~* (and prob- 
ably much closer to the entire 2.75 x 10~*) of the total flux 
in the Kepler passband can be attributed to 511B's intrinsic 
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Figure 1. Result of the simultaneous RV-photometry solution. Left: the top panel shows the filtered and phased light curve (blue) and 
the model (red), the middle panel gives a similar but expanded view, and the bottom panel includes the model residuals vs. orbital phase 
binned to f% phase. Right: the measured RVs (blue) and the model (red) are depicted vs. time (top panel), and phase (middle panel). 
The bottom panel gives the model residuals vs. orbital phase. Note phase zero is defined at mid transit, not at periastron passage. 



luminosity. A toy modej^ (using a uniform Kepler passband 
between 420 and 900 nm and blackbody spectral densities) 
gives a temperature range for 511B of Teff,BiiB = 4030°J^ 
to 4150° A" for the above contrast range. 



5 DISCUSSION 

We present the initial characterization of Kepler EB KIC 
1571511. We show that the secondary of this EB is a very 
low mass star with a mass of O.14136M0 and a radius of 
O.178317?0 - so its diameter is smaller that some planet (e.g. 
Hartman et al. 2011, Anderson et al. 2011b). For a low-mass 
object to be considered for a "benchmark status" of its class 
one would want that the object would be physically associ- 
ated with a more Sun-like star, since such stars are currently 
better understood. Better still are such binaries that are 
eclipsing, and the best constraints could come from fully- 
eclipsing and double-lined such binaries, where masses and 
radii are arrived at model-free. KIC 1571511 is almost such 
an EB - currently lacking only the secondary's RV signal. 
Indeed, if one plots the uncertainty in mass and radius of all 
well-characterized (both errors under 5%) low-mass objects 
and of 51 IB (Figure [2]|, it is easy to see that the latter oc- 
cupies a unique spot of the lowest-mass well characterized 
star just above the brown-dwarfs to stars transition (0.075 
Msun)- Importantly, the overwhelming source of error on 

see Kepler Instrument Handbook, document KSCI-19033, for 
full description. 



both the mass and radius of 511B is the error on the mass 
of the primary, 511A. This is due, on the one hand, to the 
single-lined nature of the system in the visible band, and, 
on the other hand, to the exquisite quality of the Kepler LC 
which allows for very precise determination of all LC-derived 
quantities. This, in turn, means that the continued Kepler 
observations on the target guarantee marked improvements 
in all LC-derived quantities in the future, while any obser- 
vation of the system as an SB2 (perhaps in the infrared) will 
dramatically reduce the overall error on both the mass and 
radius of 511B. This is visualized with the empty symbol of 
511B on Figure [2] which show that discounting the (model- 
ing) error on the mass of 511 A drastically reduces the errors 
on the parameters of 51 IB to potentially the best charac- 
terized and lowest-mass object - so a potential benchmark 
object indeed. 

When compared with other known low-mass stars all 
the other systems we encountered are less favorable to serv- 
ing as benchmarks of this kind: Most known late M dwarfs 
and BDs are not in eclipsing systems at all. Nearly all low- 
mass EBs, and especially those observed from the ground, 
have LCs that cannot compete with Kepler's exquisite qual- 
ity. However, there are a number other interesting objects 
already in that dataset: KOI-126 B,C (Carter et al. 2011) 
and Kepler-16B (Doyle et al. 2011) are indeed very low mass 
stars - but all are part of compact hierarchical triple sys- 
tems, making follow up and analysis more difficult. We con- 
clude that KIC 1571511B is indeed uniquely situated to be- 
come a benchmark object for the study of low-mass stars. 
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Table 3. The best-fitting model of the KIC 1571511 system and 
derived quantities. The overwhelming error source on the physical 
parameters is the error on Mi due to the SBl nature of the system 
- so the derived errors are computed twice: with Mi modeling 
errors included (marked with an asterisk) and without (statistical 
error only). 
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Figure 2. Relative precision of the mass (upper panel) and 
radius (lower panel) of all well-characterized low-mass stars 
(M* < O.55M0) with both determined to better than 5% (blue 
circles). Also noted (red square) is the location of KIC 1571511B - 
once derived only from the error on the mass of 511 A (filled sym- 
bol) and once derived only from statistical errors (empty symbol) . 
References; LHS 6343 C (Johnson et al. 2011), WASP-30b (An- 
derson et al. 2011a), KOI-126 B,C (Carter et al. 2011), CM Dra 
A & B, CU Cnc A & B, GJ 411, GJ 887 (Lopez-Morales 2007 
and references therein), GJ 3236 A (Irwin et al. 2009), LSPM 
J1112+7626 A and B (Irwin et al. 2011), Kepler-16B (Doyle et 
al. 2011) 
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